Non-linear and quantum optics of a type II OPO containing a 
birefringent element Part 2 : bright entangled beams generation 
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Abstract. We describe theoretically the quantum properties of a type-II Optical Parametric Oscillator 
containing a birefringent plate which induces a linear coupling between the orthogonally polarized signal 
and idler beams and results in phase locking between these two beams. As in a standard OPO, the signal and 
idler waves show large quantum correlations which can be measured experimentally due to the phase locking 
between the two beams. We study the influence of the waveplate on the various criteria characterizing 
quantum correlations. We show in particular that the quantum correlations can be maximized by using 
optimized quadratures. 



PACS. 42.65.Yj Optical parametric oscillators and amplifiers - 
fluctuations, quantum noise, and quantum jumps 



- 42.25.Lc Birefringence - 42.50.Lc Quantum 



1 Introduction 

The domain of quantum information with intense light 
beams is becoming more and more studied as practical 
applications are developed [1-3]. A large number of ex- 
periments rely on so-called EPR beams, that is a pair 
of fields containing quantum correlations on two differ- 
ent non commuting observables [4] . Several methods have 
been used to date. One way to generate them is to mix 
two independent squeezed beams produced by two inde- 
pendent OPOs below threshold [1,5], or to use third order 
non linearity in a fiber [3] . Another configuration consists 
in using a non degenerate OPO below threshold [6-8]. We 
propose here a new method where two bright entangled 
beams are produced directly at the output of a non de- 
generate OPO pumped above threshold. 

Optical Parametric Oscillators are very efhcient 
sources of non classical light [9-14]. In particular, it has 
been shown experimentally that the intensities of the sig- 
nal and idler beams exhibit very large quantum corre- 
lations [9,11,14]. ft is also predicted that the phase of 
the signal and idler should show quantum anticorrela- 
tions [15]. However, the measurement of these phase an- 
ticorrelations is not possible in a standard OPO in which 
frequency degenerate operation occurs only accidentally 
and is very sensitive to thermal and mechanical drifts. 
Thus, in practical experiments, the frequencies of the sig- 
nal and idler beam are different. This implies that the 
relative phase between the two fields has a rapidly vary- 
ing component. Furthermore, even if the device is actively 
stabilized on the frequency degeneracy working point [16], 
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this relative phase between signal and idler undergoes a 
diffusion process [17], similarly to the phase of a laser field. 
In order to be able to measure experimentally these phase 
anticorrelations, one can use a phase locked OPO : one 
possible technique proposed in [18] is to use a birefringent 
plate inserted inside the OPO cavity and making an an- 
gle with the axis of the non-linear crystal which leads to 
self-phase locking. This phenomenon has been studied in 
detail in [19, 20] at the classical level. It is the purpose of 
this paper to evaluate the influence of the phase locking 
mechanism on the quantum properties of the system. 

The paper is organized as follows. In sectional we give 
the equations describing the quantum fluctuations of the 
system. We study in section Ol the noise spectra of the 
amplitude quadratures difference and phase quadratures 
sum showing quantum noise reduction, i.e that these vari- 
ances are below what would be observed for independent 
coherent states. In section^] various entanglement criteria 
are discussed, namely EPR correlation properties (section 
14.1(1 , inseparability criterion (section I4.2|l and teleporta- 
tion fldelity (section |4.3|I . Finally, section [S] is devoted to 
the case of non-orthogonal quadratures which enable us 
to increase dramatically the violation of the inseparabil- 
ity criterion. 



2 Basic equations 

In order to calculate the fluctuations for the signal and 
idler fields when the OPO is pumped above threshold, 
we apply the input-output linearization technique de- 
scribed e.g. in [26]. We first determine the classical sta- 
tionary solutions of the evolution equations : these solu- 
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tions are denoted A- where i = 0, 1, 2 corresponds to the 
pump, signal and idler fields. We then linearize the evolu- 
tion equations around these stationary values by setting : 
Ai ^ Al +SAi, z = 0, 1, 2. Introducing the input fluctu- 
ations in all the modes, we obtain a set of linear coupled 
equations for the fluctuations of the different flelds. By 
solving these equations in the Fourier domain, one can 
calculate the individual noise spectra as well as the corre- 
lations between the different quadratures. 

Let us recall the system under study which has been 
discussed in detail in [20]. We consider a x'^^ medium 
placed inside a ring optical cavity doubly resonant for sig- 
nal and idler (fig.ni. Type-II phase-matching results in or- 
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Fig. 1. Ring cavity OPO with the fluctuations corresponding 
to the fluctuations coming through the input mirror (SAY^) 
and the fluctuations due to the losses (SBl"). 



thogonally polarized signal and idler fields. The reflection 
coefficients for signal and idler are assumed to be equal in 
module : rj — rexp{iQ), j = 1,2. We assume also that 
the cavity has a high finesse for these modes so that we 
can write r = 1 — k with k <C 1 ; the transmission is thus 
t « ^/2k. Finally, the cavity is characterized by the round- 
trip losses for the signal and idler waves (crystal absorp- 
tion, surface scattering), /x, and we define a generalized 



reflection coefficient : r' = r(l — /i) « 1 — {p + k) = 1 — 

Let us denote / the crystal length and L the free prop- 
agation length inside the cavity. The crystal indices of 
refraction are ni and n2 respectively for the signal and 
idler waves. The waveplate is rotated by an angle p with 
respect to the crystal neutral axes. In contrast to [20], we 
assume here that this angle is small. 

The crystal is birefringent and can be characterized by 
its mean index fi, birefringence 9 and effective nonlinear 
coupling constant g : 



ni + 712 



— (na -ni)l 



9 = Ix^ 



_ 1 (2) I fiU!QU!li02 



2c^Sanonin2 



(1) 
(2) 

(3) 



where coo/2 is the common frequency of signal and idler. 
The mean round trip phase-shift is S = ^ (ne + nl + 

L)+C where ^ is the mean mirror birefringence, C — ^^3^- 
We assume that the signal and idler modes are close to res- 
onance so that the corresponding round trip phase shifts 



s± 



+ C2 - Ci 



(4) 



where tp is the birefringence of the waveplate, are close to 
an integer multiple of 2tt : 



Sia = 27rpi,2 + ^1,2 



(5) 



where pi_2 are integers and Z\i_2 are small detunings. 

In this case, the classical evolution equations given in 
[20] become : 
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where A™ is the input pump amplitude at frequency ujq. 

The stationary solution of these equations, correspond- 
ing to the phase locked regime, can be found only within a 
given range of detunings A2) which has been detailed 
in [20]. 

The evolution equations for the fluctuations then are : 



dSAi 
dt 



SAi{-^' + iAi)+ 



g (^A^o ^SA; + A^2 ^*SAo^ + i2pe'^'^~^'^6A2+ 



d6A2 
dt 



5A, 



is) 



9 {A'^5Al 
V2^5A^ + 



\/2k(5A1" + V2Ai(5Sf 

6A2{~Pl' +iA2) + 

a['^*5Ao^ +i2pe-'^'^-'l'^SAi + 

(7) 

where (5 A™ and 6BI" correspond to the vacuum fluctua- 
tions entering the cavity due respectively to the coupling 
mirror and to the losses. From these equations, one can 
calculate the Fourier components of the fluctuations and 
thus the spectra of variances and correlations for different 
quadratures of the signal and idler flelds. 

3 Variances 

The amplitude quadrature and phase quadrature fluctua- 
tions can be written in the Fourier domain as 

6S:,{n) = 6A,{n)e~"^' + SA*{~Q)e''^' , (8) 

Sp^{n) = -I (<5i,(/?)e-*'^' - (5i*(-^2)e*'^') , (9) 

i = 1, 2 where 8Ai denotes the Fourier component and Lpi 

is the phase of the stationary solution AL'f' . 

To characterize the entanglement, one can calculate 
the normalized variances of the sum and difference of the 
amplitude and phase quadratures [6,21] 



s^{n)^\l^8~x^{n)±b~x2{n)f) 



\8p^{n)±bp2{n)\' 



(10) 
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The fact that and Sp are smaller than one is the signa- 
ture of simultaneous quantum correlations on both phase 
and amplitude quadratures. This is the case in a stan- 
dard OPO above threshold with the restriction already 
mentioned in the introduction that because of the phase 
diffusion, can not be measured easily. 

/^From the above calculated fluctuations, one can ob- 
tain the variances for any operating point. For the sake of 
simplicity, we consider here the case of detunings A — ±2p 
which correspond to the minimum threshold for which 

i — 1,2 where 



- = I and A\'^ = Vs ''^ ~ 
a is the input pump intensity normahzed to the threshold 
of the standard OPO. We now have : 



%4 + 4i72 



s+{n) = 1 



[2^ 



(11) 
(12) 



where 12 



is the frequency normalized to the 



iJT U 

cavity cut-off frequency, f2c- 

The expression for the phase sum variance is identi- 
cal to the one obtained in the case of a standard (non 
phase locked) OPO : one obtains very large phase anti- 
correlations for any waveplate angle. The expression for 
the intensity difference variance is different from the ex- 
pression obtained in the case of a standard OPO and does 
not depend on the pump power a. The intensity difference 
noise reduction is now degraded at low frequencies by the 
phase locking : inhibition of phase diffusion prevents 5*" 
from going to -|-oo as !7 goes to zero. Thus, in order to 
verify the corresponding Heisenberg inequality, Sp cannot 
go to zero but instead takes a finite value. 

The corresponding spectra are shown on fig |21 and El 
These figures show that the amplitude quadrature dif- 




Fig. 2. Intensity difference variance spectra as a function of 
the normalized frequency O for waveplate angles p — (dot- 
ted line), p = 0.01, (dashed line) and p = 0.05 (continu- 
ous line). Intracavity losses and output coupling coefficient, 
(/i, k) = (0,0.05). The standard quantum limit is shown in 
grey. These plots are independent of a. 



ference variance goes to zero as p goes to zero while the 
phase quadrature sum goes to zero close to threshold. 



Fig. 3. Phase sum variance spectra as a function of the normal- 
ized frequency O for normalized pump intensity a = 1 (dotted 
line), a = 1.1 (dashed line), a — 2 (continuous line). Intracav- 
ity losses and output coupling coefficient, {p, k) — (0,0.05). 
The standard quantum limit is shown in grey. These plots do 
not depend on the waveplate angle p. 



An important characteristic of the system is the min- 
imal value obtained for S~ as a function of f2, (S~)^- 
which is plotted in fig.^jas a function of p. This plot shows 
that one obtains a significant noise reduction as long as p 
is small compared to k. 



4 Entanglement criteria 
4.1 EPR correlations criterion 

In the Gedanken Experiment discussed in [4] , one tries to 
infer the quantum state of a particle by measurements on 
a spatially separated particle : the optical equivalent of 
that scheme is to try to infer the signal beam by mea- 
surements on the spatially separated idler beam. In this 
case, quantum correlations do not suffice to characterize 
the efficiency of the system and another criterion has to be 
considered [21]. The noise on the measured observable has 
to be taken into account : a large noise will reduce the pre- 
cision with which the unmeasured system is known. One 
is lead to introduce inference errors or conditional vari- 
ances [21]. The conditional variances for amplitude and 



(SJ™ 



Fig. 4. {S^)^^^ as a function of waveplate angle p for intra- 
cavity losses and output coupling coefficient, {p, k) = (0, 0.05). 
The standard quantum limit is shown in grey. 
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phase are given by : 



= min(|(5ii 
g 



gSx2 



mi 



Vp = min(|(5pi - 5(5p2p) = {\Spi - Cp(5p2p) 



where 



Cx — 



Cp — 



{dpi{Q)Sp2{~f2)) 



(13) 
(14) 



(15) 
(16) 



are the normalized correlations. Ref. [21] shows that one 
must have VxVp < 1 in order to have a demonstration of 
the EPR paradox. 

Even with correlations factors close but not equal to 
unity, the conditional variances can be large if the quan- 
tum noise on the individual beams is large as is the case 
for the amplitude quadrature. 

The conditional variances Vx and Vp are shown on fig. 
and |B| Due to the large excess noise on the amplitude 
quadrature at low frequencies, Vx is large for low frequen- 
cies while Vp remains low for all frequencies. This differ- 
ence between the two conditional variances is one of the 
main differences between our case and the other devices 
which have been introduced to generate bright entangled 
beams [5-8]: in these cases, one obtains theoretically sym- 
metric behaviors for amplitude and phase which is not the 
case in our system. 




Fig. 5. Conditional variance spectra for the amplitude quadra- 
ture Vx as a function of the normalized frequency O for wave 
plate angle p = (dotted line), p = 0.01, (dashed line) and 
p = 0.05 (continuous line). Intracavity losses and output cou- 
pling coefficient, {p, k) — (0,0.05). The standard quantum 
limit is shown in grey. Normalized pump intensity a = 1. 



The variation ofVxVp with J7 is plotted on fig.H This 
figure shows that EPR correlations exist in a broad fre- 
quency range even for large values of the waveplate angle. 
The behavior of the criterion is dominated by that of Vp 
which goes to zero as f2 goes to zero and is independent 
of the waveplate angle. 




Fig. 6. Conditional variance spectra for the phase quadrature 
as a function of the normalized frequency O for wave plate 
angle p = (dotted line), p — 0.01, (dashed line) and p — 0.05 
(continuous line). The first two curves cannot be distinguished 
on the figure. Intracavity losses and output coupling coeflicient, 
{p, k) — (0,0.05). The standard quantum limit is shown in 
grey. Normalized pump intensity a — 1. 
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Fig. 7. VxVp as a function of J7 for intracavity losses and 
output coupling coeflicient, {p, k) — (0, 0.05) and for wave 
plate angle p — (dotted line), p — 0.01, (dashed line) and 
p = 0.05 (continuous line) . The EPR correlations limit is shown 
in grey. Normalized pump intensity a = 1. 



4.2 Inseparability criterion 

However, this criterion is not always relevant. In the vast 
majority of quantum information protocols, what matters 
is the non-separability of the quantum state that is the 
impossibility to express it as a product of independent 
sub-states. In order to evaluate the inseparability of two 
systems, Peres proposed a necessary criterion which ap- 
plies to finite dimension systems [22] which was latter 
shown by Horodecki [23] to be sufficient for dimensions 
smaller than 2x2 - or 2x3-. In the particular case of 
Gaussian states, Duan et al. [24] and Simon [25] devel- 
oped a necessary criterion for inseparability which can be 
written as : 



Sx + < 2 



(17) 



where S~ and 5*+ are the previously introduced variances 
spectra. Let us note the entanglement produced by the 
self-phase locked OPO is not in the so-called standard 
form ( [24], equations 10-11) : in this case, the above cri- 
terion is a sufficient but not necessary condition for insep- 
arability. 
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The value of + Sp as a function of i7 is plotted 
in figlHl At low frequencies, the behavior is dominated by 
the behavior of which is large. As higher frequencies 
are considered, the figure shows that the system is indeed 
entangled over a large range of frequencies. When p <C k, 
the excess noise on the amplitude quadrature makes the 
sum larger than two meaning that the system is separable. 
The curves are plotted for a value cr = 1, that is at ex- 
act threshold. Increasing a does not increase significantly 
the variance and does not modify Sp and will thus 
move all the curves toward the inseparability limit while 
keeping the sum below two. It is the behavior of the am- 
plitude fluctuations which dominates that of the sum of 
the variances. 



Fig. 9. Fidelity, .7^ as a function of analysis frequency f2 
for intracavity losses and output coupling coefficient {fi, k) = 
(0, 0.05) and wave plate angle p = (dotted line), p — 0.01 
(dashed line), p = 0.1 (continuous line). The grey lines cor- 
respond to perfect teleportation (jT — 1), "best copy" limit 
= 2/3) and classical limit {T =1/2). 




Fig. 8. + Sp as a function of for intracavity losses and 
output coupling coefficient {/i, k) = (0, 0.05) and wave plate 
angle p = (dotted line), p — 0.01 (dashed line), p = 0.1 
(continuous line). The inseparability limit is shown in grey. 



obtains a fidelity larger that the "best copy" limit > 
2/3) over a large frequency range. 

In order to determine the optimal operating parame- 
ters, one can plot the values of waveplate angle p, normal- 
ized intensity a and analysis frequency f2 for which the 
fidelity is equal to the "best-copy" limit, that is T — |. 
The plot consists of two zones : close to threshold, i.e. 
a < 1.1, one obtains a fidelity larger than | for large val- 
ues of p and for SI < 0.5 (dark grey surface) ; for a > 1.1, 
a large fidelity is observed over a larger range of frequen- 
cies but only for very small values of p, typically p < 0.02 
(light grey surface). 



4.3 Fidelity of teleportation 

In teleportation experiments, one tries to transfer quan- 
tum states from one point to another with a unity gain. 
In this case, the successful of the teleportation is char- 
acterized by the overlap between the input and output 
states. This overlap is usually called the fidelity [1]. It can 
be shown that in order to obtain a fidelity larger than 
|, quantum correlations must be used, while a fidelity 
larger than | denotes the possession by the receiver of 
the best copy available [2] . In the particular case of gaus- 
sian states"'^, one can show that the fidelity is related to 
the intensity difference and phase sum variances defined 
in equation 110! bv the relation : 



(18) 



The fidelity as a function of the normalized frequency 
is shown in figl^jfor different angles of the waveplate. One 
sees that for small values of the waveplate angle p, one 




and only in this case 



Fig. 10. Values of the analysis frequency f2, wave plate angle p 
and normalized pump intensity a for which .7-" = | . Intracavity 
losses and output coupling coefficient (p, k) — (0, 0.1). 



One observes that in order to obtain a fidelity sig- 
nificantly larger than |, one can use very small angles, 
typically less than 1° (0.015 radians). As mentioned in 
ref. [20], this is feasible experimentally, since a small an- 
gle does not reduce the range of the crystal temperature 
and cavity length where the locking is observed but only 
its total surface : once frequency degeneracy has been ob- 
tained, existing active stabilization of the system allows 
for operation in this zone. However, working very close 
to threshold will reduce the parameter range where phase 
locking is observed and may be difficult experimentally. 
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We have shown in this section that our system pro- 
duces bright entangled beams which fulfill diflferent en- 
tanglement criteria over a large range of frequencies. At 
low frequencies, the entanglement is observed for p <C k, 
that is for a rotation of the waveplate small compared to 
the mirror transmission. 



5 Use of optimized quadratures 

The limitation of the EPR correlations zone to small angle 
is due to the fact that the intensity correlations are no 
longer perfect in a phase locked OPO as was the case in 
the standard OPO. A simple picture shows the origin of 
this phenomenon. Ai and A2 having a fixed phase relation, 
the output beam is in a fixed polarization state which 
we denote Alright ■ The orthogonal polarization mode has 
a zero mean value and is denoted Adark- One can show 
that they correspond to the polarization eigenmodes of 
the "cold" OPO, that is without the pump beam [27]. 
One can also show that their variance spectra are 

Sx{Adark) = {\5S:dark\^) = \ - 55:2?) S- 

\ 

Sp{Ab„ght) = {\5pbright?) ^ -j^{\5pi + 5p2?) ^ 

(19) 

The exact expressions of the mean fields Abright and Adark 
depend on the choice of the operating point. For the min- 
imum threshold point where we have performed our cal- 
culations, we have 



Im(E) 



Abright 
-^dark 



A1+A2 

A1-A2 
V2 



(20) 
(21) 



which are left and right circularly polarized due to the ^ 
phase shift between Ai and A2. 

As we have seen previously that S~ and Sp can be 
smaller than the corresponding standard quantum limit, 
Abright IS phasc squcczcd while Adark is squeezed on the 
orthogonal quadrature. This is shown in fig. II II which rep- 
resents the Fresnel diagram common to the fields Ai, A2, 
Abright and Adark, the phase reference being the phase of 
the field Abright ■ 

The increase in S~ at low frequencies corresponds to 
to the tilt of the corresponding squeezing ellipse (in black) 
with respect to the standard OPO situation (vertical el- 
lipse in dark grey). Due to the presence of the waveplate, 
the amplitude quadrature of Adark (measurement of the 
width along the horizontal axis) is not the optimal quadra- 
ture : it is thus natural to detect this optimal quadrature 
instead of the amplitude quadrature [28]. It is possible 
to dephase^ appropriately Ai with respect to A2 so that 
Abright and Adark are now squeezed on orthogonal quadra- 
tures. In the simple case of the minimum threshold point, 



^ for instance, by using a combination of A/2 and A/4 wave- 
plates. 




Re(E) 



Fig. 11. Fresnel picture of the signal and idler (white surfaces), 
Ahright (light grey) and Adark fields in the case of a standard 
OPO (black) and in the case of a phase locked OPO (dark 
grey). 



the new entangled are simply obtained by a linear su- 
perposition of the original fields along the neutral axis of 
the crystal. In the simple case of the minimum threshold 
point, these fields are given by : 



A' -1- 
A' 



cosi'-fU,-^sm{^\A2 



sin! ^ Ml 



I cos I I A-? 

2 ; \ 



where 6*0^4 is given by : 



^opt 



and 



1 

= — arctan 



2C— •\/ Sx (Adark) Sp {Adark) 
dark) Sp {Adark) 

{Sx (Adark) Sp {Adark)) 
\/ Sx (Adark) Sp (Adark) 



(22) 
(23) 

(24) 
(25) 



The inseparability criterion developed by Duan et al. [24] 



Sx {A'bright+) + Sp (A'dark) ^ ^ 

where 

Sx,p (Abright. dark) are defined above and A'f^right.dark 
the optimally squeezed states : 



(26) 



are 



A' 



bright. dark 



A[ ± A'2 
V2 



(27) 



Setting p and a, one can plot Sp (^A'^„g,,t^ +Sx (^^arfc) 
as a function of the frequency (figure I12|l . This figure 
shows that one can obtain a much larger violation of the 
inseparability criterion using optimized quadratures. For 
p = 0.01 rad, the value obtained is very close to the curve 
obtained for a standard OPO (where frequency non de- 
generacy and phase diffusion prevent the measurement). 
Furthermore, the divergence obtained at low frequencies is 
completely eliminated. For larger values of the waveplate 
angle, the improvement is still noticeable. 
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Fig. 12. Sp (A'ij^igftt) +Sx (A'jark) a function of analysis fre- 
quency J? for intracavity losses and output coupling coefficient 
(/i, k) = (0, 0.05) and waveplate angle p = (dotted line), 
p = 0.01 (dashed line), p = 0.05 (continuous line). Grey lines 
correspond to the curves obtained on fig.|Hlfor Sp + 



6 Conclusion 

We have shown that a type II OPO containing a birefrin- 
gent element can be a very efficient source of intense entan- 
gled beams. This entanglement has been characterized by 
the use of several criteria, namely intensity difference and 
phase sum variances, inference errors which denote the 
possibility to make an EPR-type Gedanken Experiment, 
inseparability criteria which correspond to the quality of 
the entanglement and fidelity which represent the qual- 
ity of teleportation with unity gain. These criteria have 
been shown to be fulfilled for a large range of operating 
parameters thus showing the quality of the system. The 
inseparability criterion can be largely improved by the use 
of optimized quadratures. We are currently setting up an 
experiment based on this system. 
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